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Starting from the d-wave RVB mean-field theory of Kotliar and Liu, we present a new, long- 
wavelength/low-energy exact, treatment of gauge fluctuations. The result is a theory of gapless 
fcrmion quasiparticles coupled to superconducting phase fluctuations. We will discuss the physical 
implications, and the similarity /differences with a theory of BCS pairing with phase fluctuations. 
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The high-temperature superconductors are doped 
Mott insulators. Shortly after the experimental report 
of high T c , Anderson proposed the "resonating- valence- 
bond" (RVB) idea. [Q According to this idea the root 
of high temperature superconductivity is an insulating 
spin liquid of singlet pairs. Due to Coulomb blockade 
these pairs are unable to move. Doping creates electronic 
vacancies that mobilize these pairs. Once mobile, the 
valence-bond pairs can Bose condense into a supercon- 
ducting state. 

From the beginning the Neel-ordered state found in 
the undoped cuprates presents a difficulty to the RVB 
idea. Motivated by the experimental fact that the 
Neel order is destroyed once the sample becomes metal- 
lic/superconductive, there is a general hope that suffi- 
cient doping stablizes the long-sort (doped) spin liquid. 

The notion of RVB is made very attractive due to 
the experimental data on the "underdoped cuprates" . 
For example there exists a significant temperature range 
(T* > T > T c ) where a pseudogap exists in the spin 
excitation spectrum when there is no superconductivity. 
Angle-resolved photoemission spectroscopy (ARPES) ||] 
revealed the fact that the size and the momentum de- 
pendence of the pseudogap is roughly the same as the 
d-wave gap in the superconducting state. In addition, 
ARPES studies of undoped cuprates find a peak above 
the Mott-Hubbard gap with a dispersion similar to that 
of spinons || in the RVB mean-field theories. 

On the theory side, implementation of RVB beyond 
mean-field theories is hindered by the necessity to 
treat strong gauge fluctuations. J^-||] In the absence of a 
reliable treatment of the gauge field, it is impossible to 
tell which features of the mean-field theory actually cor- 
respond to reality. In this paper we revisit the Kotliar- 
Liu d-wave RVB mean-field theory. || What we were 
able to find is an exact treatment of gauge fluctuations 
at long wavelength and low energy. Eqs. jT^Jl^ ) plus the 
discussions after it 1-6) are the main results of this paper. 

Our starting point is the following boson-fermion rep- 
resentation of the t-J model: [|| 
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Here i,j runs through the sites of a square lattice, ip., 
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the rest of the notations are standard. In the literature 
bt and /j CT are often referred to as holon and spinon an- 
nihilation operators respectively. 

First, let us briefly review the results of mean- field the- 
ory- HQ For x > 0, the mean-field solutions are char- 
acterized by the following parameters: iaoi(t) — ► /i/, 
bi(t) -> b, Xij(t) -> Xo, and A tj (t) -> A %-. Here 
rjij = +1 if i,j are are i-nearest-neighbor, and = —1 
if i,j are y-nearest-neighbor. When the temperature is 
not too high, both fif and xo are nonzero. Depending 
on the values of b and Ao the low temperature x — T 
plane is divided into four regions: i) d-wave supercon- 
ducting phase where Ao ^ 0, and 6^0, ii) the spin gap 
phase where Ao ^ 0, and 6 = 0, iii) the Fermi liquid 
phase where Ao = 0, and 6^0, and iv) the strange 
metal phase where Ao = b = 0. The boundaries between 
these four phases are Tb e (x) and T pa i r (x) at which b and 
Ao develop expectation values respectively. In mean-field 
theory both Tbe{x) and T pa i r (x) mark phase transitions. 
This statement will be made invalid by fluctuations. A 
schematic phase diagram can be found in, e.g., Ref. [ ^0| . 

According to the phase diagram jl(| T pa i r (x) > 
Tbe{x) in the underdoped regime. Throughout this 
regime there exists a temperature range Tbe(x) < T < 
Tpair (x) m which a spin gap exists without supercon- 
ductivity. In the rest of the paper, we shall concentrate 
on the d-wave superconducting phase and the spin gap 
regime mentioned above. 

At temperatures much lower than T pa i r the important 
fluctuations include: i) the gapless spinon excitations and 
ii) 
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Among these fluctuations aoi and act to enforce the 
no-double-occupancy constraint. More explicitly ani en- 
forces ftohrj + bfbi — 1, while ay ensures no net flow 
of the boson and fermion current. In the following we 
present a new approach that treats the aoi and et^ fluctu- 
ations (hence the constraint) exactly at long wavelength 
and low energy. 

First let us write down a continuum action (its micro- 
scopic cutoff length is the holon separation) to capture 
all the fluctuations discussed above: 
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In the above = \[p~b~ e = \[Pb4>b is the holon 
field, <f> S p = e l9sp is the phase factor associated with 
the spinon pairs, pb and p sp are the average den- 
sities of holon and spinon pairs, p is the density 
at half-filling, m& is the holon effective mass, is 
the physical gauge field, Cf is the Dirac quasiparti- 
cle action which depends on the pairing amplitude A , 

J f = (\zZnK^ n a^V F ^^ }a ^V F ^+^ 2a ) is the 

3-current of the spinon quasiparticles, where (n = 
1,2) are the spinon Nambu spinors associated with the 
two Dirac points respectively. We note that the spinon 
and spinon-pair parts of Eq. (||) are the action of a phase 
fluctuating BCS superconductor. |ll| 

By substituting ipb = \fpb~4>b into Eq. (0) we obtain 
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where Spb = pb~ pb- In obtaining the above we have used 
the fact that due to the lattice effect (i.e. when viewed by 
the vortices, density p of bosons can be gauged away be- 
cause it corresponds to one flux quantum per plaquette) 
ipbft&fa = -2ip sp cf>* b f cf> b . 

Next we Hubbard-Stratonavich decouple the first, the 
fourth, and the fifth terms to obtain 
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In Eq. (|) J b p = (Sp b ,jb), Jspp = (Sp sp ,jsp) where j b , 
j sp and 6 p sp are the auxiliary field introduced by the 
Hubbard-Stratonavich transformation. Physically Jbp 
and J sp p are the three-currents of the holon and spinon 
pairs. In the following, for distances greater than the 
holon separation (~ 1/ \/x) , we shall linearize Eq. (||) by 
replacing pb in the first term by pb and drop the third 
term. The resulting action read 
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In the above Kb 
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Now we can integrate out exactly. The result is the 
constraint 



Jbp + 2{Jspp + Jfp) — 0, 
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i.e. the total three current with respect to a M is zero! 
Physically this is due to the no-double-occupancy con- 
straint which requires each lattice site to be occupied by 
either a holon or a spinon. Since they both carry the 
same charge with respect to a M , the total charge is con- 
stant at any time. As the result the total three-current 
vanishes. If we define J pfl = {8p p ,'} p ) = \Jbp, Eq. (||) 
implies 

Jspp = —Jpp ~ Jfp- (9) 
Substitute Eq. (§) into Eq. (|) we obtain 
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In the above 
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Eqs( |l2| - |l4| ) is the main result of this paper. Aside from 
the renormalization factor z p and Zj and the quasiparticle 
interaction in Eq. @ the form of Eq. © agrees with 
that written down in Ref. |ll[ based on weak-coupling 
considerations. In the following we comment on several 
salient features/consequences of Eqs.(|2|-|lJ). 

1. Quasiparticle interaction: The last two terms in 
Eq. ( |l4| ) correspond to quasiparticle interaction. Since 
both of them are local and have scaling dimension 4 with 
respect to the free Dirac theory, they are irrelevant in the 
renormalization group sense. In addition both coupling 
constants „ 2 . „ and _ ± 2 , _ t stays finite as x — > 0. 

(As nOwe expect Kb ~ a;, while A sp , Mb and M sp all 
stay finite.) For these reasons we believe that one can 
safely neglect the quasiparticle interaction. 

In the superconducting phase, one can further inte- 



grate out the gaussian phase fluctuations in 



J sp- 



The 



generated spinon quasiparticle interaction is again local, 
and irrelevant. This result suggests that in the supercon- 
ducting state the spin degrees of freedom are asymptoti- 
cally described by a free Dirac theory. 

2. Quasiparticle charge and current renormali- 
zation: An important effect of the no-double-occupancy 
constraint is to introduce renormalization factors z p and 
Zj . We emphasize that when Zj ^ 1 and z p ^ 1, Eq. (|12] ) 
is not equivalent to a BCS superconductor with phase 



fluctuations. |jll|,|12[ Given the x-dependence of Kt yS p and 
Ub yS p, it is easy to show that Zj ~ x while z p ~ 1. With 
such x-dependent Zj Eq. (12) does not have the desired 

to explain the temperature 



form written down in Ref. 



derivative of superfluid density. Indeed, according to the 
present result the suppression of the superfluid density 
due to thermal quasiparticle current fluctuation will be 
proportional to x 2 T. H|iO[ It is interesting to note that 



since z p ~ 1 we expect an x-independent linear temper- 
ature correction to u p due to quasiparticle screening. 



It is worth emphasizing that in deriving Eq. ( |l2| ) 
we have not assumed holon condensation, and that 
Eqs.( p^| , p^|Jl4| ) are applicable in both p -coherent and <f> p - 
incoherent phases. 

3. Flux quantization: After integrating out a p , only 



the combination <f> p = 4>b4>sp °f the holon and the spinon- 
pair phase factors appear in Eq. (p^). We emphasize 
that the appearance of (K in the above expression is not 
due to holon pairing. This is important because in the 
presence of Coulomb interaction, and in the absence of 
screening at distances smaller than inter-holon separa- 
tion, such pairing will be highly energetically unfavor- 
able. The reason that 4>l appears in <f> p is the fact that 
the movement of a spinon pair always causes two holons 
to relocate. The combination (0* %</> p - 2^) in Eq. jl^ ) 
(note the coefficient 2 in front of ^4^) implies that in the 
superconducting state (i.e. <p p orders) magnetic flux will 
be quantized in units of hc/2e. 

4. The pseudogap phase: Under the present frame- 
work the superconducting state corresponds to the <p p - 
ordered phase, but how to view the pseudogap regime? 
According to the RVB mean-field theory, such phase 
arises from holon uncondensing. This remains qualita- 
tively true in the present results. For small x where 
the superfluid density is low, there is a wide tempera- 
ture range T c < T < T pa i r in which <j> p is thermally 
disordered while the spinons remain paired. In this tem- 
perature range there is no superconductivity but shows 
a d-wave gap for spin excitations. Clearly, it is tempt- 
ing to associate this crossover regime with the pseudogap 
regime seen experimentally. 

Due to the wide separation between T c (~ x) and T pa i r , 
there could be intermediate temperatures at which the 
thermal correlation length of 4> p is comparable with the 
microscopic cutoff ~ while is still much larger than 

the spinon pairing length ~ uf/Aq. At these tempera- 
tures the "bare" superfluid density will vanish without 
destroying the spin pseudogap. |15[ 

In the literature there exists a debate as to whether the 
pseudogap regime can be described as a phase-disordered 
superconductor. po|- p2| , |lo| Eq. (|lj) gives an affirmative 
answer to this question. However, as pointed out earlier, 
due to Zj , z p ^ 1 and the discussion below, Eq. ( |l^ ) does 
not describe an ordinary BCS superconductor with phase 
fluctuations. 

Lastly, since we have integrated out the gauge field ex- 
actly, the issue of whether the gauge fluctuations desta- 
bilize the pseudogap regime |l6| does not arise. ]i7| 

5. Spin-charge separation: Eq. ( |T^ ) exhibits spin- 
charge separation in the sense that the low-energy spinon 
excitations are largely determined by the RVB mean-field 
theory, and are little affected by the gauge fluctuations. 
On the other hand the charge response depends on quan- 
tity such as Kp which is greatly affected by the gauge fluc- 
tuations. A concrete example of this separation is that 
the vortex core size will be very different from spinon 
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pairing length scale (~ vjt/Aq). 

To 2 calculate the size of vortex core we restore the 
term in Eq. <^). Repeating the calculations be- 
tween Eq. (0) to Eq. (|lj) we obtain the following (most 
relevant) addition to Eq. (|l2|): 
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2m b p b 

Using this result we obtain a vortex core size given by 
J_ _ u sp + Au b 

e ~ 

' m b p b 

Thus as x — > the vortex core size diverge as -4=. In 
sharp contrast, the size of spinon pair (vf/Aq) does not 
depend on x. A further manifestation of the spin-charge 
separation is in the structure of vortex core described in 
the following. 

6. Vortex core: Near a vortex core p b become spa- 
tial dependent so that p&(x) — > toward the center of 
the vortex. Under such condition K b in Eq. (0) be- 



comes spatial dependent and K b (jc) 



Pb(x) 



toward 



the center of the vortex. As the result both K p and Zj 
vanishes in the vortex core. On the contrary the spinon 
excitations (which depend on the mean-field pairing am- 
plitude Ao) remain largely unaffected by the depletion of 
superfluid density. Q (To be more precise, since holons 
are depleted from the vortex core, the RVB mean-field 
pairing amplitude Ao should slightly increase toward the 
center of the vortex.) This is qualitatively different from 
the behavior found in BCS superconductors where the 
pairing amplitude (or local gap) collapse at the center of 
the vortex. This difference could be manifested in the 
absence of vortex-core induced midgap states. |HJ 

At closing it is important to point out that the treat- 
ment in this paper eventually becomes invalid as x — > 0, 
where a charge SU(2) symmetry emerges. We be- 
lieve that the enlarged gauge symmetry is responsible for 
antifcrromagnetism. Moreover, the SU(2) gauge fluctua- 
tions at length scales smaller than the holon separation 
can cause important modifications of the bare parame- 
ters (or even introducing new spinon interaction terms) 
in our theory. Finally, the cutoff length and energy scales 
in this work can be affected by the (dynamic) stripe cor- 
relation ]23j] presents in the underdoped systems. 
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